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We have developed a method for calculation of quantum fluctuation effects, in particular of 
the uncertainty zone developing at the potential curvature sign inversion, for a damped harmonic 
oscillator with arbitrary time dependence of frequency and for arbitrary temperature, within the 
Caldeira-Leggett model. The method has been applied to the calculation of the gray zone width AIx 
of Josephson-j unction balanced comparators driven by a specially designed low-impedance RSFQ 
circuit. The calculated temperature dependence of AI^ in the range 1.5 to 4.2K is in a virtually 
perfect agreement with experimental data for Nb-trilayer comparators with critical current densities 
of 1.0 and 5.5 kA/cm^, without any fitting parameters. 



The current attention to quantum information pro- 
cessing (see, e.g., the recent monograph jl)) has re- 
newed interest in fast "single-shot" quantum measure- 
ments, especially in potentially scalable solid-state sys- 
tems. Among such systems, superconductor "balanced 
comparator" , based on two similar Josephson junctions 
(Fig. la), stands apart as a very simple, scalable system 
for which quantum-limited sensitivity has already been 
demonstrated experimentally Q. 

The device is essentially a SQUID (see, e.g., |^) in 
which two similar junctions are biased in series by a 
source of Josephson phase difference ipeit), and in par- 
allel by the current Ix to be measured. Let the system 
with \(j)f,\ < TT settle in an equilibrium state (j) — 4>i, and 
then apply a rapid phase change A0e = 27r. (This can 
be readily done using the so-called RSFQ circuitry - see, 
e.g., the recent review Q.) As a result, the system be- 
comes statically unstable and the Josephson phase cj) has 
to switch to one of adjacent stable states, depending on 
the sign of I^- (For junctions with substantial damping, 
the choice is limited by two states closest to the initial 
value of 0: (/) f ~ (pi ± tt) . 

This process may be readily understood using the 
" magnetic language" : the driver circuit providing the 
pulse A(j)e — 2TT in fact injects a single flux quantum 
into a superconducting loop formed by its output stage 
and the comparator (Fig. la). Since the loop is low- 
inductive (non-quantizing), the flux quantum has to drop 
out across one of the comparator junctions, depending on 
the sign of Ix- This transient process produces a large 
(discrete) output signal, the so-called SFQ pulse V{t) 
with JV{t)dt = = h/2e across the corresponding 
junction. Such a pulse may be readily picked up and 
registered by relatively crude devices so that the ac- 
curacy of the Ix sign measurement is defined entirely by 
the comparator. 

In the absence of fluctuations, the boundary between 
the two possible outcomes would be infinitely sharp - see 
the dashed line in Fig. Ic; however, fiuctuations create a 
finite " gray zone" of Ix where probability p of switching 
to a certain finite state changes gradually from to 1 - 
see the solid line in Fig. Ic. The gray zone width, which 
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FIG. 1: (a) The balanced comparator, (b) its potential en- 
ergy profile U{(j)) for Ix = 0.01/c and two values of the ex- 
ternal phase <j>e (/i = cos(<^e/2)), and (c) probability p of 
comparator switching into state 0/ = 0i — vr as a function of 
the measured current Ix (schematically). 



is traditionally defined as 
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characterizes the accuracy of the single-shot measure- 
ment. This width, including its temperature dependence, 
has been measured in several special experiments with 
externally-shunted Josephson junctions P,0]. However, 
theoretically it has been only calculated 19,0 for a spe- 
cial function 0e (t) enabling an analytical solution of the 
problem, but rather different from that used in experi- 
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ments. Thus, the comparison of theory and experiments 
was not completely conclusive. The goal of this work 
has been to develop a general method of calculation of 
A/^ for an arbitrary waveform (jj^^ (t) and to compare the 
results with experimental data |5| . 

The potential energy of the balanced comparator (Fig. 
la) may be presented in the form 
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where Ej = hlc/'2e is the Josephson coupling energy 
scale, and Ic is the critical current of a single junction. (In 
the simplest case, neither (t) nor depend on the state 
of the comparator.) Notice that the part of "washboard" 
potential profile, contributed by the Josephson junctions, 
changes sign when 4>e{t) is increased beyond tt. This is 
exactly the reason of the system switching to one of the 
newly stable states 0/ ~ 0i ± tt (Fig. lb). Now let \Ix\ 
and the natural current scales of thermal and quantum 
fluctuations. It — {2e/h)T and Iq = {2e/h)huj„ = 2eujp, 
respectively (where ujp is the plasma frequency ) , all be 
much smaller than I^, and the potential inversion time 
be of the order of, or shorter than the characteristic time 
of system dynamics. Then the choice of the final state is 
determined by the system evolution close to the point — 
0. In order to describe this evolution, we may keep only 
two leading, linear and quadratic, terms in the Taylor 
expansion of the potential energy (0) near this point: 
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This means that the state choice problem in the orig- 
inal, nonlinear system is reduced to that of a damped 
time-dependent harmonic oscillator with frequency coit) 
defined as uj'^{t) = n{t)ujp, where fi is switched rapidly 
from a positive initial value fii (in experiments 
close to 1) to a negative final value ^/ = — /i^. 

The probability of switching to a final state with 0/ < 
4>i may be found as 
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where p{4>,(j)',t) is the system's density matrix traced 
over the degrees of freedom of the environment, and 
4'max{t) is the coordinate of the maximum of potential 
(^ after inversion. Converting to coordinates rj = + (j)' 
and ^ = (/) — (/)', we can express p(0, 0, t) = {l/2)p{r], 0, t) 
via the system propagator J{r],^,t\r]i,^i,0): 

-t-oo 

p(77,0,i)= J{r],0,t\r],,^,,0)p{r],,^,,0)dr]id^i. (5) 



To find the propagator, we may use the Caldeira- 
Leggett approach |^ with the linear distribution of the 
environment oscillators, which gives a quantitatively 
correct description of systems with externally shunted 
Josephson junctions. According to this theory. 
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X J dr J ds^(r)^(s) cos[w(r — s)]. 




Here, C = M{4))'^/2 - Mfj.{t)uj'^4>'^ /2 is the Lagrangian of 
the mechanical oscillator equivalent to our system, with 
mass M = 2Ep/ujp, while 7 — ujp/2uJc is its damping 
parameter, where lOc = {2e/K)IcR, and R is the shunt- 
ing resistance Parameter is the cutoff frequency 
of the environment oscillators. To evaluate the path in- 
tegral (||), we represent coordinates 77,^ as a sum of the 
path parts 7?(T),^(r) minimizing the action 5, and small 
fluctuations ry, ^. The path parts satisfy the following 
equations M: 
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It is convenient to present solutions of these equations 
as follows: 



r){T,t) = r]iai{T,t)+T]a2{T,t) + {I^/Ic)aiT,t), (11) 
ar,t) = Mr,t)+Cb2{T,t), (12) 

where functions ai.2 and 61.2 as functions of r obey the 
uniform versions of equations (^, ( [lO| ) with the following 
boundary conditions: ai(0) = &i(0) = 1, ai{t) = bi{t) — 
0, 02(0) = 62(0) = 0, a2it) = b2{t) = I, while a is the 
solution to Eq. (^ with the unit right-hand part and 
boundary conditions a(0) — a{t) — 0. 

Plugging all these expressions into Eq. (H), and car- 
rying out a lengthy but straightforward (Gaussian) inte- 
gration, we get 
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where F'^it) is a normalization factor, and the prime rep- 
resents differentiation over r. 

These formulas present the generalization of Eq. (6.26) 
of Ref. Q to the case of arbitrary time dependence of the 
oscillator potential curviture /i(i). Equations (^), ( p^ ) 
show that if the initial density matrix is Gaussian (as it 
is, e.g., for a system in thermal equilibrium), the final 
matrix is also Gaussian, with the average phase ((/>) and 
variance (0^) determined by parameters Ki,N,Qi and 
C. (Other parameters affect only the final phase velocity 
distribution, which is not important for our particular 
problem.) 

Using the definition ([l]) , the gray zone width may now 
be calculated as 
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This is our central result. For the particular case of the 
RSFQ driver circuit used in experiments |j , the func- 
tion iJi{t) has been calculated numerically from the cir- 
cuit schematics, using the PSCAN software package 
Since functions a(T, t) and 6(t, t) are exponential near the 
boundary points t — Q and r ^ t, the standard "shoot- 
ing" methods for the numerical calculation of these func- 
tions would be unstable. Because of this we have used the 
relaxation method jlOj. Upon the calculation of a(T, t) 
and 6(t, t), parameters Ki and Qi were obtained by the 
standard numerical integration using the trapezoidal ap- 
proximation . C was calculated using 3D Monte Carlo 
integration where an integration by parts helps control 
the discontinuity at a; — > 0. Due to the shape of the func- 



tion bi , a combination of stratified and importance sam- 
pling greatly increases the convergence time, so a variant 
of the VEGAS algorithm |jll| was used. The bath oscil- 
lator cutoff frequency f2 was taken large enough (SOwp) 
to avoid any effect on the calculation results. 

Figure ^ shows /i and the essential parameters of the 
Gaussian distribution as functions of t. One can see that 
if the interval [0, t] includes the time point ij„„ of the po- 
tential curvature sign inversion, with both ti„„ and t—tinv 
much longer than the oscillator's reciprocal bandwidth 
Au^^ « max[l/a;c, t^c/2cjp], then C > 1, and Qi » Ki, 
so that the final density matrix and switching probability 
p do not depend on the initial state of the system. In this 
limit, Eq. (fsl) takes a very simple form: 
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Figure |^ shows the resulting temperature dependence 
of the gray zone width for several values of the iner- 
tia parameter (normalized junction capacitance) Pc = 
((jjc/wp)^. At high temperatures, AI^ grows as T^^-^ due 
to thermal fluctuations, while at T ^ it saturates due 
to quantum fluctuations. Note also that the dependence 
of AIx on Pc is different for high and low temperatures: 
if thermal fluctuations dominate, the gray zone width 
depends on Pc only weakly, saturating at comparable 
values at both /?c — > and Pc 00. However, in the 
quantum fluctuation range {T —> 0), AIx grows as /Jc^^ 
at high damping (/3c — > 0) and saturates in the oppo- 
site limit of low damping. All these dependences may be 
qualitatively understood from the following simple con- 
sideration: AIx crudely equals to the signal current that 
creates the phase shift (j) = Ix/2Ic equal to the r.m.s. 
value of phase noise in thermal equilibrium. The lat- 
ter value may be estimated assuming that an equivalent 
current noise source H with equilibrium spectral den- 
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FIG. 2: Calculated function fi{t) and parameters of the final 
phase distribution for the RSFQ drivers used in experiments 
^ for f3c — 1. The time scale uj~^ is close to 1.1 ps for the 
critical current density jc ~ 1 kA/cm'^ |^ and 0.47 ps for jc 
= 5.5 kA/cm^ 11. 
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FIG. 3: Temperature dependence of A/a; for shown 
in Fig. 2, calculated for several values of fi^- Dashed lines 
represent the thermal limit. 



sity S'/(cj) = (4/i?)(?iw/2) coth(?icj/2r) acts on a time- 
independent linear oscillator within the bandwidth Acj 
defined above. 

Figure ^ shows the comparison of our results 
with experimental data for comparators based on 
niobium-trilayer (Nb/A10a;/Nb) Josephson junctions 
with /c |t=4.2K = 145 [lA^ Pc — 1, for two values of 
the critical current density: jc — 1 kA/cm^ |Q and 5.5 
kA/cm^ Q. One can see that besides the deviation of 
the two lowest-T points in experiments Q , which was ap- 
parently caused by sample self-heating, the theory gives 
a virtually perfect description of experimental results, 
without any fitting parameters. (The possibility of a 
substantial external noise contribution to A/^; in experi- 
ments [|[ H] has been ruled out by special control experi- 
ments using similar comparators, fabricated on the same 
chip, but driven by "softer" waveforms.) 
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FIG. 4: Temperature dependence of A/x. Points show exper- 
imental data from Refs. ^. The dashed lines are results 
of calculation taking into account the Ambegoakar-Baratoff 
temperature dependence of Ic- The solid lime is the theory 
for an instantaneous change of /x from -1-1 to -1. 



To summarize, we have developed a method of anal- 
ysis of quantum fluctuations at the inversion of the po- 
tential curvature sign of a damped harmonic oscillator. 
When applied to Josephson junction comparator, these 
results may be used for numerical caculation of the gray 
zone width A/^,. Such calculation for the Nb-trilayer 
comparators |^ gave a nearly perfect agreement with 
experimental data. Our result may be also generalized 
to the case of a finite inductive impedance of the source 
of the signal Ix, which is typical for Josephson junction 
systems, e.g., magnetic flux qubits Indeed, the 

impedance may be described by connecting the source 
inductance L in parallel with the source of Ix in Fig. 
la. An elementary calculation shows that this leads to 
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This means that if the inductance is not too low, A > 1/2, 
the input SFQ pulse A0e = 27r develops an instability of 
phase 4> just as was described above, and our theory gives 
a ready recipe for the calculation of /S.Ix and hence of the 
signal flux resolution /S.^x — L/S.Ix- 

However, in order to reduce dephasing, flux qubits typ- 
ically require unshunted Josephson junctions. This is 
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why a natural next task would be a calculation of 
for the case when damping is dominated by quasiparticle 
tunneling in unshunted junctions. For this, the Caldeira- 
Legget action should be replaced with one found by 
Ambegaokar et al. p^ . 
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